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I. INTRODUCTION 



Modelling the scattering problem of nucleons by nuclei using a complex-valued potential 
leads to a basic physical model applicable in the region of low energies 0, 0| . This approach 
is called the optical model and it is often used for the analysis of experimental data. The 
parameters of the complex- valued potential are fitted by solving the associated Schrodinger 
equation and then attempting to obtain the best agreement with the experimental scattering 
cross-section. This method requires laborious numerical calculations 

In this paper we develop another model, called the zero-range potential with inner struc- 
ture, which serves as a generalization of the Fermi zero-range potential model describing 
low-energy neutron scattering by nuclei . The Fermi zero-range potential is generalization 
by adding a non-trivial inner Hamiltonian and appropriate inner space 0] supplied with an 
indefinite metric. This zero-range model with inner structure effectively takes into account 
multiparticle interactions of the scattering particles which is important for such problems as 
scattering of electrons by atoms and scattering of nucleons by nuclei. The zero-range model 
with inner structure is a solvable model which is an important feature of our approach. This 
is the main advantage of our model over the optical model. Moreover we do not introduce 
a non-Hermitian Hamiltonian as in the optical model, hence the S-scattering matrix of our 
model is unitary and has the correct placement of poles and zeros in the complex plane of 
the wave number k. 

It is known also [5] that the descriptions of the resonance processes and collisions with 
redistribution are the most difficult problems in the optical model. In the present work 
we consider the resonance scattering in the s-channel where all results are given in explicit 
analytical form. In this section we examine the main works starting from the fundamental 
paper which laid the foundation for the development of the zero-range model. It is useful 
to study the Fermi model, assuming that the nucleus potential is 3-D finite square well: 

v { r )= r v °' iir - r ° (1) 

0, if r > r 

where Vq > 0. Let us consider the bound state of the particle in this well with the energy 
E = — f^ft 2 < 0, where fi = 7 ^ 1 1 ^ 2 is the reduced mass of scattering particles. Then the 
solution of the Schrodinger equation decreasing as r — > oo is: 

C 

ip{r) — — exp(— nr). (2) 

Defining the function x(r) = rip(r) and using Eq. (J2J) we find the condition on the boundary 
of the 3-D square well : 

1 d , A 

-X{r) = -k. (3) 



X(r)dr' , , r 

Let us consider the special case, i.e. the weakly bound state: \E\ « V or equivalently 
ktq « 1. This holds, for example, in the case of deuteron and negatively charged ions. If 
in this case k, or the energy of the bound state E < 0, is fixed and the inequalities m" << 1 
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and krQ << 1 hold, where k is the wave number, then the boundary condition (J3J) may be 
reduced to the condition : 

T^l-X(r)) =-k. (4) 
.X{r)dr J r=Q 

This boundary condition was introduced by Fermi in the scattering problem of neutrons 
by nuclei and means that the wave function ip(r) is singular at r = and in the vicinity of 
the zero point r = has the form: 

V>(r) « C (~ - kV (5) 

Introducing the Fermi s-wave scattering length a = k~ x one can find that if A; — ^ then the 
total cross-section will be o = Ana 2 . The boundary condition (@J) determines the so called 
the Fermi zero-range potential model [3]. It is important that the boundary condition (0J) 
holds in the low-energy limit k — > for a deep well with arbitrary form when the bound 
state with the energy E = -g K 2 < is a weakly bound state. 

We note that, in the general case, the Laplace operator is symmetric on the scattering 
waves of the Fermi zero-range potential. Really, if the asymptotic of the scattering waves at 
the origin r = is: 

^(r) = A + b, + (1), (6) 

then the corresponding boundary form for general smooth elements *&i(r) possessing the 
above asymptotic at the origin r = can be written 

J ik = (- A %, - (%, - A v&fc) = AiB k - A k Bi, (7) 

where 

= y ^ t (r)^ fc (r)d 3 r. 

Due to Eq. (JSj) the values Ai and 5j satisfy to the boundary condition 

kA { + ATiBi = 0, (8) 

hence the boundary form (jUJ) is zero: Ji k = 0. 

In the case of square well when the conditions K,r << 1 and kr << 1 hold one can make 
a small correction to the Fermi result : 

- = «-y«, 0) 
where a -1 = — (A;cot5(/c))/ c= o and the function kcotS(k) has the form : 

kcot5(k) = -- + jk 2 . (10) 

Here S(k) is the scattering phase in the s-channell (1 = 0). Eqs. Q and (jlOj) are also correct 
in the low-energy limit when the well has an arbitrary shape but the top bound state with 
the energy E = —j^k 2 < is the weakly bound state. In this case the parameter tq is the 
effective radius of the theory and is only in an approximate sense the well radius Ji|. 
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A single standard Fermi zero-range potential defines an operator with only trivial spectral 
structure (one negative eigenvalue). Still the question remained until recently: is it possible 
to extend the construction of Fermi to include operators with rich spectrum? A version of 
this question appeared in the paper [(J by Wiegner, where it was noticed that for compactly 
supported potentials V(r) (V(r) = for r > r ) the scattering waves may be calculated for 
r > r if we use an appropriate energy-dependent boundary condition on the sphere r = r : 

9 X (r,k)+a(k)x(r,k)] = 0, (11) 



dr 

' I —I u 

where x{ r i k) = ri/)(r,k). Unfortunately this condition contains the full scattering data, 
which correspond to compactely supported potentials. Really, the solution of the Schrodinger 
equation leads to the s-scattering waves at r > r for compactely supported potentials in 
the form : 

C 

ijj(r,k) = — sin(fcr + 8(k)), (12) 
r 

where S(k) is the scattering phase in the s-channell (1 = 0). Combining Eqs. (|llj) and (|12|) 
we find : 

a(k) = -kcot(kr + 5(k)). (13) 

Hence to calculate a(k) we should know 5(k) which is equivalent to the solution of the 
Schrodinger equation. We note that the suggestion of Wigner inspired numerous attempts 
to construct the energy-dependent potentials and appropriate boundary conditions, see also 

HI- 

The authors of the paper |9( suggested a new version of the zero-range potential in terms of 
von-Neumann Operator-Extension techniques. All books and papers on zero-range solvable 
models in quantum mechanics, see for instance 0,1 

lllj . are based on this important paper 

Thereby Gelfand attracted [12] attention of mathematicians to the necessity of reconsid- 
ering the main results of von-Neumann's Operator-Extension theory for general Hermitian 
operators in terms of symplectic Hermitian forms similar to the boundary form of a differ- 
ential operator. The above form (J7J), for instance in Ri, is: 

\x)ip k {x) + i)i{x)i)l(x))<lx = ($(z)V>k(a:) - ^i(x)ip' k (x)) x=0 . (14) 

Realization of this plan would imply some sort of "integration by parts" for abstract Her- 
mitian operators, see for instance section 2 below. The first step in the development of 
symplectic techniques in Operator-Extension theory was done in 0] with the intention of 
creating a quantum mechanical solvable model with a reasonably rich and algebraically 
computable discrete spectrum and resonances. This program was developed in a series of 
publications, see for instance [ill EH EH| • 

In |4| the "zero-range potentials with inner structure" were suggested as an application 
of the classical Krein formula [16. [T7j | to "spectral modelling" of quantum systems. This 
prospect was anticipated by M. Krein himself long before [l8|. In particular the operator 
extension approach permits one to construct simple solvable models of quantum systems 
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with resonance properties without solving sophisticated boundary problems, but operating 
with finite matrices (see also and more references in fill]). 

Unfortunately zero-range solvable models, in particular models with inner structure, con- 
tain a large number of parameters which have no physical interpretation and hence are not 
fitting parameters. We propose a solution of this problem in this paper which has application 
for low-energy resonance scattering of neutron by nuclei. The generalization of the Fermi 
zero-range potential is obtained by a symplectic version of the operator-extention technique 
and assuming that the zero-range potential has inner structure which arises by adding some 
"inner" Hamiltonian and "inner" space E with a special indefinite metric. Though our treat- 
ment is based on Operator-Extention techniques we are able to define generalized zero-range 
potentials by special boundary conditions at the origin r = 0, which is a generalization of 
the Fermi condition (£Q) : 



This condition can be found formally if we suppose that the scattering wave has the form 
(I12|) at the origin r = — which is not trivial because, for conventional potentials, it is not 
generally true. It is important to understand that this boundary condition does not define 
the particular zero-range potential but it defines a whole class of zero-range potentials. In 
section 5 we show that the zero-range potentials derived in our paper do in fact satisfy the 
boundary condition fT5|) . 

We note that the boundary condition (fTKj) follows from the expantion of the wave function 
ijj(r, k) at the origin r = as: 



Another, most important, point of our treatment is the principle of analyticity of the 
function: 



Thus we require that this function is an entire function of the complex variable k which 
permits us to evaluate all the extra parameters of the zero-range model with inner structure. 

In the next section we develop symplectic techniques for the general extension procedure 
of symmetric operators in a Pontryagin space (this may be thought of as an abstract version 
of integration by parts). In the third section we calculate the S-scattering matrix based 
on a variant of the Krein formula. Then in section 4 we study a special class of zero- 
range models of scattering systems for which the scattering matrices fulfil the analyticity 
condition formulated for the function kcot5(k). It appears that for the models of this class, 
if the scattering length, effective radius and the eigenvalues of the inner Hamiltonian are 
fixed, then all other parameters of the model, in particular the boundary parameters and 
the deficiency vector in the finite- dimensional case dimE = N may be calculated directly. 
Moreover, the type of the indefinite metric in the inner space E is also pre-defined by the 
above condition of analyticity. 

Then in section 5 we apply the results to a typical problem of resonance s-scattering. In 




(15) 




(16) 




(17) 



n=0 
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Appendix A we give a straightforward derivation of the basic Krein formula. One usually 
assumes that the simplest zero-range potential is the singular potential proportional to the 
Dirac 5-function, which is often used in various physics problems. Consequently, we provide 
additional on this case in Appendix B. 



II. SYMPLECTIC VERSION OF THE OPERATOR-EXTENSION TECHNIQUES 

At first we define the dimensionless coordinates X{ in i? 3 space by the equation x = k r as- 
suming that ko is some characteristic wave number. We also assume that the wave functions 
■u(x) are dimensionless then the Schrodinger equation for scattering problem can be written: 

(A + A - V"(x)) u(x) = 0, (18) 

where A = E/E is a dimensionless energy with E = and V"(x) is the dimensionless 
potential 

\/(x) = £ -V(r), E = ^. (19) 

Here \i = mim 2 /(mi + m 2 ) is the reduced mass of the scattering particles. Below we will 
use conventional notations 

k 2 

x:=x = k r, A := -2, (20) 
then the equation for the Greens function of free evolution is: 

- (A + X)G x (x,x') = 5 {3 \x-x'), (21) 
and the appropriate Greens function for outgoing waves has the form: 

G x (x, x') = -(A + X)- l 5^\x - x 1 ) = X \ . (22) 

Att\x — x I 

We develop here the symplectic technique for the zero-range potential of a pair of quantum 
particles in three-dimensional space, assuming that in the course of interaction the pair may 
enter an "inner" space E supplied with an indefinite metric. This interection is determined 
by an inner Hamiltonian B, which is a Hermitian operator with respect to the properly 
defined dot-product [. , .] in E. The two-body Hamiltonian in the outer space, with respect 
to coordinates defined in terms of the center of mass of the pair, may be given by a three- 
dimensional self-adjoint "non -perturbed" Schrodinger operator in the "outer" space L 2 (i? 3 ) 
with the conventional dot-product (u,av) = a{u,v) = (au,v): 

lu = -Au, (23) 

defined on the Sobolev class W%{Rz) of all square-integrable functions u possessing square 
integrable second-order derivatives. In the course of our analysis the spectral parameter 
may take negative and complex values. The operator I = — A is defined with respect to 
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the dimensionless space variable x. Other variables x , y will be used for vectors in the inner 
space E of the model, see below. 

Following we restrict the operator Z to the domain consisting of all smooth functions 
vanishing near x = 0. The closure Zo of the restricted operator is defined on the domain 
Do which consists of all VFf -functions vanishing just at the origin x — 0. This operator is 
symmetric and has deficiency indices (1,1): 



(Z - U)D = L 2 (R 3 ) e Ni, (Z + U)D = L 2 (R 3 ) e AL i; 

dim Nj = dim N_i = 1. The adjoint operator Z^ is defined, due to the von- Neumann theorem, 
see 22j, by the same differential expression (J23j) on the domain Dq = D + Ni + iV_j, where 
the one-dimensional subspaces iVj,iV_j are spanned by the corresponding non-perturbed 
Green functions G\(x, 0) at \m\f~X > and X = ±i : 

Ni = {G-i(*,0)}, N^ = {G 4 (*,0)}, 



-i\x\ p iV~i\x\ lA-jT 

G-,(x, 0) = ^-j-p, Gl (x, 0) = j^ = Y-jG-iix, 0), (24) 

where we define the branch of the square root by the condition Im\/A > 0. 

The above description of the domain of the adjoint operator based on the von-Neumann 
formulae was used in J^, |2^. The other description of the domain of the adjoint operator 
as the set of singular elements described above (JHJ) was probably invented by Fermi [3] and 
is now commonly used in the physical literature, see for instance Each element u 

from the domain of the adjoint operator is characterized by the asymptotic boundary values 
A, B for |x| -> 0: 

A u 

u(x) = —— + B u + o(l). (25) 

It is well known that these representations are equivalent. Nevertheless we have a good 
reason, see Lemma 2.2, to derive this fact from the decomposition of elements of the defect 
Ni + iV„j with respect to the symplectic basis 

W + (x) = i [G^(x, 0) + Gi(x, 0)] = J^G-iix, 0), 

W-[x) = i [G_i(x,0) - GfaO)] = - — T — G?_i(x,0). 

We will present the von-Neumann formula for elements of the domain of the adjoint operator 
in terms of a decomposition with respect to the basis W± : 



u 



uo + 77+W+ + 7£W_, (26) 



where Uq is an element from the domain of the closure Zo of the restricted operator and 
?7± are complex coefficients. The deficiency elements G±i(*,0) are the eigenvectors of the 
adjoint operator with eigenvalues =F« respectively. Then the elements W± are transformed 
by the adjoint operator as: 

l+W + = W-, Z+W_ = -W+. (27) 
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Lemma II. 1 The boundary form 

(l£u, v) - (u, l£v) := Ji(u,v), 

of the adjoint operator 1$ depends only on the terms u d , v d of the elements u, v in the defect. 
It is an Hermitian symplectic form in the variables r]± and may be presented as 

Ji{u,v) = Ji(u d ,v d ) = (l£u d , v d ) - (u d , l£v d ) = 

-L= (rf + ft - rj-fjl) • (28) 

The symplectic variables rj± are connected with the asymptotic boundary values A, B via the 
transformation: 

x 10 

-1 -AnV2 

Then the Hermitian symplectic form is given by Ji(u,v) = ( K B U A V — A U B V> ). 







\ n- 1 





(29) 



Remark The above representation of the "outer" boundary form in terms of the sym- 
plectic variables A, B is usually obtained, see for instance 0, E3, Eif via straightforward 
integration by parts. We use the construction suggested in the above lemma to introduce 
new symplectic variables which are useful in the symplectic version of operator-extension 
theory for general "abstract" operators, see below. 

Proof of the formula (|2*Hj) may be obtained via direct application of the adjoint operator 
to the above version (|2*Bj) of the von- Neumann decomposition of elements from Dq\ 

l+u = l u + V u + W„- V u _W + , 

in the boundary form of the adjoint operator: 

(IqU, v) - (u, l£v) := Ji{u,v). (30) 

One can easily see that the above "outer" boundary form depends only on the parts of 
the elements u, v in the sum of the deficiency subspaces (on the "defect"), Ud = r]+W + + 
ri u _W-, v d = r]'" + W+ + rf_W-: 

Ji(u,v) = Ji(u d ,v d ). 
Now the boundary form may be calculated (using (|27jl ) : 

(l^Ud, v d ) - (u d , l+v d ) = (vlC - V u _f)l) [ \G^,0)\ 2 d 3 Z. 

The integral is evaluated as l/(47rv2) giving us the desired formula (|28J) . The second formula 
is based on the calculation of the asymptotic behavior of the non-perturbed Greens function 
as x — > : 

w + (x) =ik-^ + ° (l) ' w - {x)= -^ +0(1) - 



8 



Hence for x — > we have the asymptotic behavior: 

rj + W + (x) + v-W-{x) = rri- ! T ± 7 = t + o(l). 

Att\x\ Any/2 

The second announced formula follows immediately from the above asymptotics. 
□ 

Self-adjoint extensions of Zo were obtained in by requiring that the asymptotic 

boundary values A, B satisfy Fermi boundary conditions B = 'yA parametrized by a real 
7 = 7. In [4] an orthogonal sum of two operators Iq © B with an abstract operator B 
was restricted and then extended based on an abstract version of the above symplectic 
"integration by parts" (|2K|) in a Hilbert space E with positive metric, B : E —>■ E. We will 
now develop the symplectic technique for a g-Hermitian operator B in a Hilbert Space space 
with an indefinite metric. 

Assume that P± are two complementary orthogonal projections in a finite-dimensional 
Hilbert space E, dimE = N with conventional dot-product (*, *) and consider the operator 
g = P + — P_ defining an indefinite dot product (g-dot product): 

(gx, y) = [x,y] = (P + x,P + y) - (P_5,P_g). 

Here [ax, y] = a[x,y] = [x,ay]. We denote by E g the space E supplied with the g-dot 
product. It is a Pontryagin space with indefinite metric g if dim P_ is finite, see for instance 
[24L Ei|. We develop below a symplectic version of the operator-extension technique in the 
Pontryagin space which differs from the above techniques in L 2 (Rz) by minor details only. 

A bounded or densely defined closed operator B is called ^-symmetric (symmetric with 
respect to a g-dot product) if [Bx, y] = [x, By] at least on the domain of B in E. This 
condition is equivalent to the condition gB C B + g, where B + is the adjoint operator with 
respect to the conventional dot-product (*,*) in E. Without loss of generality we may 
assume that the operator B — il is invertible and the inverse is bounded. 

Eventually we will apply the symplectic scheme to bounded operators which are sym- 
metric with respect to the conventional dot-product in E and commute with g in the finite- 
dimensional space E (dimP = N). We assume that the spectrum of the operator B consists 
of a finite number of simple real eigenvalues X s : 

(B - X s ) e s = 0. 



One may choose a normalized generating vector |36[ e G P, [e, e] = 1 such that the 
vectors e and e' = j|rffe are connected by the above g-isometry ||^je and form a linearly- 
independent pair (Be 7^ 0). We define the domain Pb of the restricted operator as Pb = 
[B — il]~ l {E Q g {e}} and set B = B| D . Here the orthogonal difference with respect to 
the g-dot product is denoted by Q g . The vector e plays, in this construction, the role of the 
deficiency vector of the restricted operator at the spectral point i: 

[(B-i7)D Bo> e]=0, 

and the vector e' ([e',e'] = [e,e]) plays a role of the deficiency vector at the spectral point 
— i: 

[(B+iI)D Ba ,e>] = 0. 
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The sum M of two one-dimensional subspaces Mi, M^i spanned by the vectors e, e' respec- 
tively is called the defect of the operator B . At any other two complex conjugated regular 
points A, A of the operator B the deficiency subspaces N x , N x are one- dimensional. The 
corresponding deficiency vectors are calculated as 

B + il B + il 

e x = — r-e, e- x 



B — A/ B-AJ ' 

and the defect M\ + M\ is two-dimensional. Let us define the g-adjoint operator B^ by the 
formula 

[B x,y] = [x,B$y\, 

on the elements y for which the g-dot product in the left hand side may be continued onto 
the whole space E as an "anti-linear" functional [i?} of x. For the densely-defined operator 
B and e chosen from the complement of the domain this condition implies 

(B+ + iJ)e = 0, (B+-iI)e' = 0. (31) 

For a bounded operator B we just define the adjoint operator on the defect by the above 
formula (}3*Tj) . Then for any complex value of the spectral parameter the deficiency vector is 
the eigenvector of the adjoint operator at the adjoint spectral point: 

(B+ - A/) e A = 0, (B+ - A/) e- x = 0. 

This is similar to the above fact for Laplacian : 

IqG\ — AGa, 

on the space of square integrable functions L 2 (R-z). 

An analog of the von-Neumann representation of the domain of the adjoint operator 
remains true in a Ponryagin space, see |24(. So, for elements from the domain of the adjoint 
operator we have the representation x = x + ae + /3e', x G -Do with 

Bq (x + ae + (3e') = Bx Q — iae + i/3e'. 

Note that in the case when the original operator B commutes with g and is self-adjoint 
with respect to the conventional dot-product, the operator is isometric both in the 
conventional and g-metric. 

Consider a new basis for the defect M = Mj + M_j of the operator B , which is similar 
to the above basis consisting of linear combinations of Greens-functions: 

e + e' B e-e' -I 

w + = = — -e, W- = = — -e. 

+ 2 B-il ' 2i B-il 

We see from fjBTjl that 

BqW + = BqW- = —w + . 

In the symplectic version of operator-extension theory we use the new basis to represent 
elements from the domain of the adjoint operator via new dimensionless symplectic variables 
t^± which play a role analogous to the above pair r]±: 

x = x + £lw + + 
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Lemma II. 2 The boundary form of the adjoint operator may be given by an Hermitian 
symplectic form in terms of the variables £±: 

K(x, y) = [B+x, y] - [x, B + y] = CC - C ■ (32) 

Proof exactly follows the pattern of the previous lemma concerning the representation of 
the boundary form in the "outer" space L2CR3). Really, the value of the boundary form 
depends only on the parts of the vectors x, y in the defect, 

K(x,y) = 

[b+(£>+ + cw-), (&w + + e.w-)\ - [(£>+ + cw-), b+(c> + + ew-)] = 
k^«;_ - cw + ), ie + w + + e.w-)] - [(£>+ + cw-), - e.w + )] = 

C+e([u-, w-) + K, u;+]) - Ce(K, u;_] + K, u;+]) = ^ - 
since the coefficients ([w-,w + ] — [w + ,wJ\) in front of Q£± are equal to zero. Thus we have 

B —I 

w< = — — e, W- 



B-il ' B-i/ 

which implies ([iy_,w_] + [tu + , w+]) = 1 and ([w-.,w + ] — iy_]) = and the announced 
abstract integration by parts formula (|32p. 

□ 

Consider now the orthogonal sum of operators /©B in the Pontryagin space L 2 {Rz) ®E g 
with elements Z7 = (m, x), m G L,2(R), x G Restricting both operators as above we 
obtain the symmetric operator Iq © Bo with defect iV © M and deficiency index (2, 2). We 
define the adjoint operator as an orthogonal sum of adjoints Iq © Bq and calculate the 
boundary form of it as a sum of boundary forms of the outer and inner components: 

J(U, V) = J(u, v) + K(x, y) = (B U A V - A U B V ) + - ^M), 

where we choose the symplectic variables A, B for the "outer" boundary form. The bound- 
ary form J~(U, V) of the orthogonal sum of operators Iq © B^ in the Pontryagin space 
L 2 (B g Eg is a symplectic Hermitian form. Hermitian extensions of the operator l © B (or 
Hermitian restrictions of the operator Iq © Bq") may be constructed via imposing boundary 
conditions on symplectic variables A, B,£± such that the total boundary form vanishes on 
the corresponding planes. For instance the condition defined by a Hermitian matrix is 



r = r 



7oo 7oi 
7io 7n 



B / 7oo 7oi M 

-C- / \ 7io 7io / I £+ 



(33) 



The planes in the four- dimensional space of complex symplectic variables A, B,£± where 
these boundary conditions are fulfilled and hence the boundary form vanishes, are called 



11 



Lagrangian planes of the symplectic form J~(U, V). We assume here that the elements 7^ 
have numerical values. 

All Lagrangian planes of the above form may be constructed either with a help of a 
Hermitian matrix T, as above, or obtained from the constructed plane by a proper j7-unitary 
transformation, for instance: 

A — ► A', B — ► B', 

which leaves the form J{U, V) unchanged. One may show that all Lagrangian planes of the 
form J~(U, V) may be obtained via the composition of Hermitian transformations (|33|) and 
the fundamental ^/-unitary transformation A' = A, B' = B, = = — 

The matrix elements 7^ are called "the boundary parameters ". They do not have 
any straightforward physical interpretation. One of the most important problems of the 
theory of solvable models based on operator-extension methods is the problem of fitting 
these parameters in a physically meaningful way. We solve this problem for the general case 
in section 5. 



III. KREIN FORMULA AND THE S-SCATTERING MATRIX 

The expression for the non-stationary S-scattering matrix obtained as a product of wave- 
operators may be derived |2jj from the corresponding Krein formula llTT l . This formula 
gives a description of all resolvents of an extended operator in terms of the boundary condi- 
tions (|33|) and some functional parameter — the so called Krein's Q-function, see [22J • In this 
paper we focus on the derivation of an expression for the stationary S-scattering matrix, just 
fitting the properly constructed ansatz for scattered waves to the above boundary condition 
(|33|) . In the course of the calculation of the coefficient in front of the divergent wave in the 
ansatz (the scattering amplitude) we will use the Q-function of the inner operator B^J". This 
parameter appears in the course of the solution of the homogeneous equation for the adjoint 
operator as a map of symplectic variables. The Krein's Q-function is an abstract analog of 
the Weyl-Titchmarsh m-function [jfij ]. 

Lemma III.l The Q-function for the operator B has the form: 

r I+AB . 

[e 'B3A/ e] ' 

where the symplectic coordinates £± of the solution e\ of the adjoint homogeneous equation 
are related by 

c r /+AB u 

Proof. Consider the ^-orthogonal projections P g = e][e and I — P g in the ^-metric onto 
the one-dimensional deficiency subspace N{ and onto the complementary subspace EQ g Ni, 
respectively. Then the solution of the adjoint homogeneous equation e\ may be written as: 

B + il B / . J + AB . 

e\ = e = e H e e, e + 

A B-AJ B-il B-il L 'B-AJ J 
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B^ 1 < 34 » 
The last term in the right-hand side is an element uq from the domain of the restricted 
operator Bo, but the first two terms belong to the defect M so the whole linear combination 
may be written as 

B + ^/ m 
___e = w+- Q{X)w^ + u , 

where Q(A) = [e, gge]. 
□ 

Remark The Krein's Q-function for a Hermitian operator in a space with a positive 
metric g > is a dimensionless Nevanlinna-class function with positive imaginary part in 

N 

the upper half-plane ImA > . If B is a finite diagonal matrix diag{A s } s=1 , N = dimE, 
and the metric tensor g is trivial g — I, then 



s 1 2 



L A s -A 



where the A s are non-negative eigenvalues of B and \e s \ are the squares of the moduli of the 
components e s of the deficiency vector e with respect to the (standard) basis of eigenvectors 
of B in E. 

Now consider an indefinite metric tensor g defined by a diagonal matrix, for instance 
{g ss } = ±1, s — 1, 2, . . . N, a real positive diagonal matrix B = {X s } , s — 1, 2, . . . dinrE 
and a ^-normalized deficiency vector e ([e, e] = 1) with non-zero components e s with respect 
to the same standard basis. Then the corresponding Q-function will have a form 

s=l s s=l s 

with real coefficients P s = ^ s |e s | 2 and non- negative eigenvalues A s . The Q-function (|53j) has 
poles of first order at the eigenvalues \ s of B. 

The previous lemma permits us to solve the adjoint non-homogeneous equation and derive 
the Krein formula for the resolvent of the self-adjoint (or at least self-adjoint in the g-metric) 
extension A of the symmetric operator Iq © Bo which is defined by the boundary condition 
(I33|) (see the simple derivation of the Krein formula for the operator B in Appendix A). The 
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S-scattering matrix may be derived from it in a rather standard way, see for instance 
We concentrate now on the direct derivation of an expression for the scattering amplitudes 
and the construction of special solutions — scattered waves — of the homogeneous equation 

= © B+) * = A* , 

which satisfy the boundary condition (|33|) . The scattered waves serve as eigenfunctions of 
the self-adjoint extension A of the restricted operator l © B and may be found in the form 
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where the ansatz for the "outer" component ip of the scattered wave in R3 is a sum of 
incoming and outgoing waves: 

p iV~\\x\ 

M*, v) = + T ^ ^)lTTTi ( 36 ) 

4tt\x\ 

with the amplitude T(u, y/X) in front of the outgoing wave. We note that this dimensionless 
amplitude T is proportional to the standard "physical" amplitude f : T — 4irk f. The 
component ip E of the scattered wave in the inner space E is just proportional to the limit 
value of the solution e\ of the adjoint homogeneous equation on the real axis of the spectral 
parameter: 

ip E = T E ex = T E [w+ - Q(\)w- + u ] . 
One may show, see for instance floL that the boundary conditions (|33|) formulated for 



elements of the domain of the extension of the operator 1$ © Bo are valid for the correspond- 
ing scattered waves. The symplectic variables in the above ansatz for components of the 
scattered wave in the outer and the inner spaces are: 

B=(l + i^T^j , A = T,^ = -QT E , = T E . (37) 

This gives the following equation for the amplitudes T and T E : 



(38) 



l + = / 700 701 W T 

QT E ) \j w 7n ) \T E 

Solving this equation we obtain the following expressions for these amplitudes. 

Theorem III.l The amplitudes T, and T E for real positive values of the spectral variable 
A are equal to 

where 



, A, -A 



ss ■ 



Proof may be obtained from the previous discussion. 
Remark Note that the amplitude 

T{VX) = — ' 



I701I 2 :VX 



Til— Q 47r 



contains essential information on the spectral properties of the extension A. For example, 
the solution of the corresponding homogeneous equation 

AS? = A*, 
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yields the equation 




/ \ 7io 7u 

which has a solution \1/ with T E = ^Jq\ X ) ^ anc ^ om y ^ the dispersion equation is fulfilled: 

^n/A _ _ |7oi| 2 
4vr " 7o ° 7n - Q(X)' 

The negative roots A s < of the above equation on the spectral sheet ImV^X > of the 
variable y/X correspond to the square- integrable first component of the solution hence 
they are the eigenvalues of the extended operator A. 



IV. A SPECIAL CLASS OF ZERO-RANGE POTENTIALS WITH INNER 
STRUCTURE 

In this section we study a special class of zero-range models for which the function F(k) = 
kcotS(k) satisfies an analyticity condition with respect to the wave number k ([k] = cm -1 ). 
We will show that for models from this class the boundary parameters and the deficiency 
vector in the finite-dimensional case dimE = N may be calculated directly. Moreover, the 
type of indefinite metric in the inner space E is also pre-determined by the above condition 
of analyticity. 

We will also use the dimensional physical scattering amplitude f(k) with orbital angular 
momentum I = 0, see |27j ]. connected to the above dimensionless amplitude T(y/X) in the 
"outer" space. Using Eqs.(|35|) and we find this s-scattering amplitude f(k) as: 

The scattering matrix S(k) = exp(2i5(k)) where S(k) is the scattering phase in the s-channel 
has the form: 

hence the scattering amplitude f(k) is: 

= totfw-tt " (42) 
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Combining Eqs. (j40j) . (|41j) and ()42|) one may find the S-scattering matrix in the form: 



_ , , 2ik F(k)+ik 

S(k) - 1 + Mobw -J^ j] -ik-F(k)-ik' (43) 

with 

F(k)=Arrk [ loo - l701 ''j - (44) 
7n - Q{a) 

We may rewrite the Q-function calculated in lemma 3.1 in terms of the wave number k, 
assuming that all eigenvalues X s = (k s /k ) 2 of the inner Hamiltonian B are positive. Then 
it takes the following form 

1 JL k 4 4- k 2 k 2 

QW = T^1^W P " (45) 

ft s=l K s ft 

where the weights P s = g ss \e s \ 2 are not necessary positive! Note that we have introduced the 
resonance values k s > of the wave number k through the equation X s = (k s /k Q ) 2 , hence 
[k s ] = cm' 1 . Using Eqs. (J44|) and (J45|) we obtain a similar representation for the function 
F(k) 

F(k) = kcot5(k) = 47T7oofco ^° N l1 ^ +kik T ^- 

7n ^ Z^s=i k'i-k 2 ^ s 



We consider a special type of potential with inner structure for which the function F(k) 
defined by the above formula is an entire function of the variable k. 

The behavior of the function F(k) defines the scattering process and can be represented 
by a Taylor expansion at the origin k = 0: 

^ oo 

F(k) = -- + r p 2 + ... = Y,9nk 2n , (46) 

n=0 

where a is the scattering length and ro is the effective radius. 

We may also impose on the parameters of our model the requirement that the scattering 
matrix satisfy the following "physically reasonable" behaviour: S(k) — > 1 when k 2 — > oo. 
This condition is fulfilled if and only if the denominator 7 n — Q(X) tends to zero as 1/A 
when A = (k/k ) 2 — > oo which yields: 



N '„2 



7n 



We will show that for the introduced class of zero-range potentials the S-scattering matrix 
satisfies this "physically reasonable" condition for large energy automatically. 



Lemma IV. 1 The function 

1 7oi I 



F(k) = kcot5(k) = Ank 



2 



7oo 



1 r Af fcg+fc|fc 2 p 
111 ~kj l^s=l k'i-k 2 rs 



(48) 
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is an entire function of k (a polynomial of degree N in k 2 ) 



k 2N 



F(k) « ^TrlToxi^-l)^-^, k^oo, 

with non-zero normalization constant A ; if and only if the values of all weights P s are 
defined by the formula: 

k 2N+2 A 

Here the product Ilt(t^s) ^ s s P rea d over all indices t, 1 < t < N, except one, t ^ s. This 
analyticity requirement automatically yields the boundary parameter 711 as 

J^k 2 
»=i ft o 

Proof. We may consider the denominator D of the F-function as a function of the variable 
A = k 2 /k^ with parameters A s = k 2 jk\. 

N 1+A 2 

D(X) = 7ll - Q(X) = Tii + X) X * P * ~ X \^f Ps > (50) 

s=l s=l s 

where iV = dimE. The analyticity of the fraction on the whole plane A means that there 
are no zeroes of the denominator in any compact domain of the complex A plane. Hence 
D(X), being a rational function with iV simple poles must have only one zero of multiplicity 
iV at infinity and hence both 1/D(X) and F(X) are polynomials of degree N: 

D(A)« A(_i)^ A ^00, 

where we introduce the sign-factor (— 1)^ to give the final formula a simple form. Then the 
asymptotic of the polynomial at infinity is defined by the highest order coefficient : 

A^O, F(k)^-4nk \ l01 \ 2 (-lf T . 

The constant A plays an essential role below. We will call it the normalization constant. 

On the other hand, expanding the fractions in the above representation (|5()jl of D in 
Laurent series we obtain D(X) in another form: 

N N N 

D{X) = 7 H + X * P ° + \ E( X + A ') P * + Y2 i> + A ') A ^ 

s=l s=l s=l 

N 

+ ■ ■ ■ + E(! + ^s-'Ps + o(X~ N ) « (-l)A, A - 00. (51) 

s=l 
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This means that all coefficients in front of the powers A 1 (0 < I < N) are equal to zero and 

— N 

the coefficient in front of A is equal to the normalization constant. Thus the variables 
q s = (1 + A^)P S (s = 1,2, ... ,N) fulfill the linear system obtained by comparison of the 
coefficients in front of the powers -rr (/ = 1,2, ...,N) on the left and right hand sides 
of the last equation ()51|) . This system has a positive Vandermonde determinant W = 
W(\\, A2, A3, . . . Ajv) if the parameters A s are arranged in order of increasing Ai < A2 < 
A3 < . . . Aat. The solution of it may be found as a ratio of the determinant and minors 
W s = W(\ u A 2 , . . . , A a _i, A s+1 ...X N )>0: 

(1 + A^ = (-1) S ^A= A, S = l,2,...iV. (52) 

The announced result is obtained from the last formula by inserting wave numbers k s , k 
instead of X s , A. 
□ 

Corollary 1 We may assume (see theorem 4.1), that the metric tensor is defined as 
g ss = (— l) s sgnA. Then the squares of the components |e s | 2 of the deficiency vector e are 
given by 

(53) 

Remark. The basic formula for the Q-function was derived under the assumption that 
the deficiency vector e is non-degenerate with respect to the indefinite metric form and 
[e, e] = 1. In the "generic case" 

trh u '"if a i 740 - (54) 

so we may choose the normalization constant as 

A 1 1 



,^12 _ 1 |A| 




k\kf N+1) 


1 + As Ht(t^s) 


As — A< 


(kg + kf) rit(^s) 


^s H \ 



1 tDl ^iv fc 2(JV+1) (-i) 



s=1 (fef+fcg)n t (^ S )i*?-fc2i 

and use the normalized deficiency vector in the formula for the Q-function. We assume that 
the non- degeneracy condition |5^| ) is fulfilled. 

Theorem IV. 1 The necessary condition for the analyticity of the function F(k) = kcotS(k) 
on a whole complex plane ofk for zero-range potentials with inner structure is the alternation 
of the diagonal components of the metric tensor g: 

9ss = (-l) s signA, (55) 

where the order of the components (s = 1, 2, . . .) is defined by the ordering of the parameters 
k s : k\ < k\ < — The only parameters which remain free in the model are the scattering 
length, the effective radius, and the eigenvalues of the inner Hamiltonian which commutes 
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with the metric tensor. The boundary parameters 700, |7oi| ore defined as functions of the 
scattering length a, effective radius r®, and k 2 (s = 1, 2, . . . , N.): 



~ 1 - / N \ N 



2 



In particular, if the effective radius is positive, r > 0, then A > and g ss = (— l) s . All other 
essential parameters of the model: the moduli of the components of the deficiency vector and 
the the inner boundary parameter 711 are uniquely defined from the condition of analyticity 
of the F -function on the whole k plane. 

Proof. Using Eq. (J49|) and our condition k\ < k\ < . . . one see that 

signP, = (— l) s signA. (56) 

Combining this equation with our definition P s = g ss \e s \ 2 we find the necessary condition for 
the metric tensor g given by Eq. (j55|) . The explicit form of the function F(k) follows from 
Eqs.flHJ) and flUJ): 

N 

F{k) = kcot5{k) = e - 7 Y[{k 2 s - k 2 ), (57) 

s=l 

where we have introduced the notations 

a 1 47r l7oi| 2 , KQ s 

£ = 47T7 00 fco, 7 = , 2N-1 A " ^ 58 ) 
fc A 

In particular, using (|57p. (|58|) and definition (J46|) we find the scattering length a and effective 
radius r : 

N \ _1 N 

^+7n^ 2 ,r =2 7 ^n fc *- ( 59 ) 

s=l / s=l t(p£s) 

This completes the theorem. 
□ 

Corollary 2 It is remarkable that, using (|43|) and (J57j) . we may write the scattering 
matrix 5* in terms of the parameters k s , e and 7 as: 



2ik 

e- ik -i]\ N s=l {k 2 s - k 2 ) 



S ( k ) = 1 + ~ ,.0 (60) 



The total scattering cross-section is cr(k) = 47r|/(A;)| 2 or can be written in explicit form 
(see next setion) using equations (f4"2"j) and (|H7j) . This expression for the S-scattering ma- 
trix describes the resonance scattering of particles with resonances defined by the spectral 
properties of the inner Hamiltonian (again see the next section). 
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V. RESONANCE S-SCATTERING ON ZERO-RANGE POTENTIALS WITH IN- 



NER STRUCTURE 

In this last section we consider the main physical properties of the resonance scattering 
problem on zero-range potentials with inner structure. This problem has application mainly 
to low-energy neutron scattering by nuclei. 

First we will show that zero-range potentials with inner structure, as defined above, fulfil 
our condition for generalized zero-range potentials (|15|). ie. the special boundary condition 
at the origin r = 0. Due to the fact that the condition (fTo^) follows from the expansion (fTH|) 
we will consider Eq. (j!6|) which can be rewritten in dimensionless form ()25j) as: 

u(x,k) = ^ + B(k) + o(l), (61) 

where 

A(k) = 4nk Csin5(k), B(k) = kCcos5(k). (62) 
Hence the special boundary condition (fi3j) is equivalent to the equation: 

4nB(k) k 



A(k) k 



cot5(k). (63) 



On the other hand, as follows from Eqs. ()37|) and ()39|) . the left side of Eq. ()63|) for zero-range 
potentials with inner structure is: 



AnB(k) 47T • [\ a- ( h 
A{k) = A{k) ' 



(64) 



Using Eq. we find that the right side of Eq. (J6l|) is k^Fik) = (k/ko)cot8(k), hence 
equations (jBHJ) and (fTHjl are satisfied for zero-range potentials with inner structure. Thus 
our definition (fT3j) is correct for the zero-range potentials derived in this paper. 

Let us show that our resonance scattering model depends only the scattering length a, 
effective radius r , and the spectrum k s (s = 1,2, .. . , N) of the inner Hamiltonian. It is 
convenient to define the typical wave-number k$ and a dimensionless parameter a by: 

' 47r|7oi| 2 27V-1 , R _v 

a := (47T7oo) , (65) 



47T7ooao A 

where ao is a new parameter with dimension ([a ] = cm). Then the function F(k) = 
k cot 5(k) for our zero-range potential is: 

1 N 

F(A0 = ----nW*?-«3* 2 )- ( 66 ) 



Defining the scattering length a from the equation F(0) = —a , we obtain: 

a 



1 + a nf =1 alkl 



(67) 
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Hence the parameter ao is the scattering length for the zero-range potential without inner 
structure, i.e. when a = 0. Evidently the above equation (jHTj) is the renormalization 
equation since it expresses the renormalized scattering length a via the non-renormalized 
scattering length ao, taking into account resonance scattering. 

Thus the function F(k) given by Eq. (JSfij) depends on N + 2 parameters: ao, a and 
k s , s = 1, 2, ... N. The effective radius may be found from as 



r 



N N 



2 aa oYi n a o^- ( 6g ) 

n=l s(s^n) 

Consequently, the non-renormalized scattering length a , the dimensionless parameter a and 
the spectrum k s (s = 1, 2, . . . , N) of the inner Hamiltonian define the scattering length a and 
effective radius r of the model. Vice versa, we may choose the renormalized scattering length 
a, the effective radius r , and the spectrum k s (s = 1, 2, . . . , N) of the inner Hamiltonian as 
independent parameters, and define ao and a by Eqs.(|H7jl and (}68|) but the former option is 
more useful and will be used below. 

The total cross-section for spherically-symmetric scattering is calculated as 

which implies, due to (|66j). the explicit formula: 

47ra 2 

agtf + 2a rr =1 (a§*2 -^fc 2 ) + « 2 flti (« - ag* 2 ) 



<r(fc) = ^ ^ (70) 



The latter equation exhibits clear resonance properties: for k = k s the products in the 
denominator vanish and the whole cross-section is reduced to 

which coincides with the Wiegner formula cr(fc) = j^^p for = fc s . 
Consider now the analytic structure of the scattering matrix (pfHj) 

6 ^~F(ifc)-iife 

in the complex k plane. Then the equation for the poles of the scattering matrix is: 

F(k)-ik = 0. 

One may see that all solutions of this equation may be written as k = in, where k are the 
zeroes of another polynomial with real coefficients: 



N 

a 

8=1 



II ( a o^ + o k2 ) - a K + 1 = 0. (71) 
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In the simplest case N = one may see from (|66|). (|67|). and (|71|) : 



1 nm 1 

ao = a, k = — , F(k) = — = —k, 
a a 

which is in full agreement with the Fermi model, see [2( and the Introduction above. From 
Eq. (J70|) it follows that the integral cross-section in this case is given by the classical Wiegner 
formula: a{k) = frpp > where a = a . 
For N — 1 one may find from (|71)1: 

?V - « + i = (72) 



with 



1 1 , 72 r o 

— = h aaofci, — = cta . 

a a 2 



Hence Eq. (j§f|j) in this case yields: 



F(fc) = -i + ^ 2 . 



These results and Eq. ()70l) for AT = 1 coincide with well known results for low-energy 
resonance scattering, see [2^ and Eqs. © and (jlUjl. 

In the case N = 2 one may show from equations (|7T|). (J57j) . and (|55j) that: 

r ° 4 , r 2 ,1 n fr 7Q \ 

2{kl + klf ~" + a =0 ' (?3) 

with 

1 1 9 . O 7 9 S / 7 9 7 2\ 

- = h aa k l k 2 , — = «a„ I k 1 + k 2 ) . 

a ao 2 

In this case the total cross-section is given by Eq. (fTOj) for iV = 2. In the general case N > 2 
it follows from Eq. (J71j) that some of the poles of the S-scattering matrix are situated on 
the real axis k (imaginary axis k = in = p + iq) and the other sit at complex-conjugated 
points k s = q s — ip s , R s = q s + ip s , i.e. symmetric with respect to the real axis in fc-plane, 
because all coefficients of the algebraic equation (fTTj) are real. The zeroes of the S-scattering 
matrix are defined by the adjoint equation F(k) +ik — 0, hence if a pole of the S-scattering 
matrix is situated at the point k s = p s + iq s then the appropriate zero is situated at the 
complex-conjugated point k = p s — iq s . Note that these are universal properties of the 
S-scattering matrix and are a consequence of unitarity and causality, but we have derived 
them using only Eq. (|71J1 . The imaginary poles k s = iq s may lie on both the positive or 
negative imaginary semi-axes (q s > or q s < 0), but all other poles (when p s ^ 0) should 
lie in the lower half-plane. This property is connected with general spectral properties of 
selfadjoint operators in spaces with indefinite metric, see for instance [2J, |2Jj |2j| . Thus the 
properties of the function F{k) given by Eq. (jHHjl guarantee the correct analytic structure 
of the S-scattering matrix in the complex plane of the wave-number k. 

For our model the poles k n = p n + iq n of the S-scattering matrix may be classified 
according to their positions using the following scheme: 
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1. If p n = then for q n > the pole corresponds to a bound state with the energy 
E n = — ^q\ and for q n < the pole corresponds to a virtual state with energy E n = —^q\- 

2. Poles with p n > 0, q n < correspond to metastable states with complex energy 
E n = E' n - i^f, where: 



and the rate of decay of the metastable state is w n = Here E' n > if p n > \q n \. 

3. Poles with p n < 0, g n < correspond to resonance trapping. These poles and the 
corresponding poles of the metastable states are symmetrically situated with respect to the 
imaginary axis. 

No other poles are present in our zero-range model since all poles may be found as zeroes 
of a polynomial with real coefficients. 

In conclusion we point out the main problems which should be examined in the future 
in connection with the model of zero-range potentials with inner structure. Besides pure 
resonance scattering which we have considered in this paper, scattering theory of nucleons 
by nuclei must include: 

Potential scattering connected with consideration of the interaction around the nucleus 
surface and with the Coulomb interaction in the case of proton scattering by nuclei. The 
spin-orbit interaction describing polarization effects. The generalization of the model for 
scattering with an arbitrary quantum number I. 
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VII. APPENDIX A: DERIVATION OF THE KREIN FORMULA WITH SYM- 
PLECTIC COORDINATES 

In this Appendix, using the notations B Q = B and B = B, we give a simple derivation 
of the classical Krein formula for the ^-symmetric operator Bo in the Pontryagin space E. 
This derivation follows the pattern of [l4| where similar formulae for operators in spaces 
with positive metric is presented and may be easily modified for extensions of the operator 
Z © B in Pontryagin space. 



K 



2/i 



[{Pnf ~ (qn) 2 } , r 



n 
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Consider a g-hermitian operator B in the Pontryagin space E with an indefinite dot- 
product [.,.] = (<?.,.). Assume that the operators B ± il are invertible, and £ G E is 
a finite-dimensional subspace which will play the role of the deficiency subspace Mj for 
the reduced operator B defined on D = j^jjE, Mi = {x : [(B — il) Do,x] =0}. The 
deficiency subspace at the point —i may be obtained from £ via a g-unitary transformation 
M_j = -§zff£- If the operator B is densely defined then the subspaces form a positive 
angle. If the operator B is not densely defined, (for instance, if B is bounded) then we 
assume that the angle between M±j is positive. In both cases one may define an adjoint of 
the formally self- adjoint operator, Bq, at least on the defect M = Mj + M_,. In fact the 
extension of the reduced operator B must be constructed on the defect and hence we may 
derive the Krein formula below only on the defect. 

The subspace £ is mapped into the defect by the transformations 

e — > -^i- = ^ n e:=w +1 

„ B+il „ 

e— B _ t] e _-y 

e > 2i = B-ii e - = W ~> 

where 

Bq W + = W-, Bq = —W + . 

Then the solution u = uq + -g^y£+ — -g^y£_of the non-homogeneous adjoint equation 

B + u = Buq - — Ly - yy^-y£- = Am + / 
B — il B — il 

may be found from the condition: 

We require that the symplectic coordinates £± satisfy the following boundary condition in 
terms of the g-Hermitian operator T: 

e- + = o, 

and find an explicit formula for the resolvent of the of the corresponding extension Br- 
Assuming that the operator B — XI is invertible we may multiply the above equation by 
-§ZX, t nus obtaining 

(u -n B-U I + \B C 
{B - tI)Uo= B^xi f+ B^xi^ + ^ 

Substituting here £_ from the above boundary condition and projecting orthogonally (in 
the Pontryagin space) onto the subspace £ we obtain an important connection between the 
symplectic coordinate £ + of the solution u and /: 

i+ = T-P £ ^P £ B^XI f - (74) 
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This connection permits us to calculate the resolvent of the extension: 

u= [Br-XI] 1 f = 

I + XB I c I rc I f B I 

B^Xl B^^ + ~B^^ + + ¥^Jl } + B^u^W^U L ^ 

1 :/+ ,„ [tf n .» t + = 



B-XI J (B-XI)(B-iI) 
I B + iI I B-il 

" TD \ T n I+XB n D \tJ' V^/ 



B-Xr B-XI T-P £ ^P £ B-XI 



The last expression gives the Krein formula for the resolvent of the extension B T of the 
reduced operator. 



VIII. APPENDIX B: SINGULAR 5-POTENTIALS 

We consider here potentials which are often used as an approximation of, for example, a 
hard sphere potential or a 3-dimensional deep square well potential when the atom-atom 
interactions are effectively weak and dominated by elastic s-wave scattering: 

V(r-r') = U Q 5(r-r'), (76) 

where 

TT 2 ™& 1 1 1 

U = , - = — + —. 77 

fi fi rrii f^2 

Here \i is the reduced mass (/i = m/2 when the scattering particles are identical) and a is 
the s-wave scattering length. The sign of the parameter a (a > or a < for repulsive 
and attractive interactions respectively) depends sensitively on the precise details of the 
interatomic potential. This potential is considered a good approximation when three-body 
processes can be neglected and hence we need take only binary collisions into account. Note 
that expression (|77jl for Uq is found by the Born approximation in the case of low-energy 
(k — > 0) scattering (though this approximation is is not valid here, see below). Moreover one 
can suppose that the potential (j76j) yields the total cross-section as a = Ana 2 in the limit 
k — > when the scattering particles are different and o = 8ira 2 for identical particles (Bose 
statistics). Evidently this potential is pointlike and hence, in some sense, one may consider 
it as a simple example of a zero-range potential. For Uq < this potential represents a deep 
well, see the Introduction of this paper. Nevertheless we show below that the potential (fTBj) 
is not the correct approximation for atom-atom collisions in both cases a > and a < and 
hence cannot represent some kind of zero-range potential. Actually this fact is well known. 
However, we examine it in detail as it demonstrates some non-trivial aspects of the theory 
of zero-range potentials and because this 5-singular potential is still very popular (3o| . 
We may consider the Dirac <5-function in (|75|) as a 5-sequence 5 £ (r — r') given by: 

SAr-r') { lf |r ~ r ' - £T ° . (78) 

0, if |r — r'| > er 



25 



Then the potential (|75jl is V(r) = Vq if r < etq, and V(0) = if r > etq, where 



1/0 = (79) 



In case 1 when C/q > (a > 0) and /c — > the total cross-section is [271 ]: 



/ tanh( g g r ) _ ^ 



2 



^^H^r^-V ■ (80) 

where q = y/2[iV h~ 1 ~ e~ 3 / 2 . 

For e — > we have cr ~ e 2 and hence the total cross-section ex for the singular potential 
(J7HJ) is zero which means that no scattering occurs in this case! 

We note that the conventional derivation of the Gross-Pitaevskii equation is based on 
the potential (|76|). Fortunately this equation may be found 31] without using the 5- 
approximation, see also |32^ . 

In case 2 when Uq < (a < 0) the potential (|7fi|) represents an infinitely deep well: 
V(r) = — | Vo| if r < er , and V(0) = if r > er . Then the total cross-section for k — > is: 

.2 2 Z' tang(ego^o) 



' = ^ r H~ " 7 ' (81) 

with g — a/2/x| Vol^ -1 ~ e~ 3 / 2 . The limit e — > does not exist in this equation which means 
that the potential (|7fj)l is not correct in the case Uq < (a < 0) at all. 

Thus the singular potential (|75j) is not a proper approximation for scattering processes in 
both cases 1 and 2. Actually the correct approximation of the hard sphere potential (case 
1) is given in the papers [i^ and [Hif, see also ji^f, and the correct approximation for the 
deep well (case 2) was developed by Fermi see also the Introduction. 
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